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Abstract Phylogenetic relationships among 37 Japanese skipper species were analyzed by compar- 
ing 787 bp of the mitochondrial gene encoding NADH dehydrogenase subunit 5 (ND5) and 1,110 
bp of the cytochrome oxidase subunit 1 gene (COT). Phylogenetic trees of all 37 species of skippers 
representing four subfamilies, including two pyrgine and seven hesperiine genus groups sensu 
Evans, were constructed using maximum parsimony (MP) and maximum likelihood (ML) methods. 
Our results mostly support the phylogenetic inferences of genus groups proposed by Evans based 
on morphological data. 
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Introduction 


The family Hesperiidae consists of some 3,500 species of skippers classified into more than 
500 genera, and is the only representative of the superfamily Hesperioidea (Ackery er al., 
1998). The family is considered a monophyletic group characterized by several autapomor- 
phies, such as antennae at least twice the width of the scape apart, compound eyes with a 
complete marginal ring of reduced ommatidial facets, and the conspicuously constricted 
cervical area in larvae (Ackery et al., 1998). 


A comprehensive higher classification of the Hesperiidae was first proposed by Watson 
(1893), and his system was followed by Evans (1937, 1949, 1951-1955) with some 
modifications. One of his important modifications was the proposal of the genus group al- 
most equivalent to the tribe as a taxonomic category. This introduction of the genus group 
made hesperiid classification more systematic. The subfamily classification has also been 
changed since Evans; six subfamilies are now recognized, the Coeliadinae, Pyrrhopiginae, 
Pyrginae, Heteropterinae, Trapezitinae and Hesperiinae (Ackery et al., 1998). The 
Heteropterinae were separated from the Hesperiinae, and the Megathyminae were included 
in the Hesperiinae. While molecular systematics has been occasionally applied to the 
Hesperiidae (Weller et al., 1996), the current classification of the family is based on mor- 
phological characters, including the antennae (length, shape of apical portion, nudum seg- 
ment), labial palpi (shape of palpomeres, vestiture), legs (epiphysis, tibial spinules), and 
male genitalia. Wing venation is one of the most informative morphological characters for 
lepidopteran classification and has also been used for this group. However, since wing ve- 
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nation is extremely uniform in the Hesperiidae, only the length of the discoidal cell or the 
relative position of vein origins on the cell are useful characters (e. g. Evans, 1949). 
Moreover, phylogenetic analyses of intrafamilial relationships based on morphological 
characters are rare (de Jong, 1983, 1986; de Jong & Treadaway, 1993, Chiba, 1995), and 
only a limited number of taxa has been analyzed. 


Given the current state of hesperiid classification, it is desirable to use phylogenetic analy- 
ses based on both morphological characters and molecular data to test the monophyly of 
current subfamilies and genus groups proposed by Evans. The objective of this study was 
to clarify the monophyly of some genus groups proposed by Evans, based on 37 Japanese 
hesperiid species representing nine genus groups and four out of six known subfamilies. 
Although the materials are from a limited faunal area and the molecular data were derived 
only from the ND5 and COI genes of mitochondrial DNA (mtDNA), the results are never- 
theless informative for phylogenetic systematics of Hesperiidae. 


Materials and methods 


Materials 


All 37 hesperiid species, including the recently introduced Thymelicus lineola 
(Ochsenheimer) (Naganuma, 1999), were collected in Japan (Table 1). For those species 
protected nationally or locally, official collecting permits were obtained. For 
Carterocephalus palaemon (Pallas), we collected the subspecies akaishianus Fujioka, which 
is unprotected in Yamanashi Prefecture. Our preliminary global analyses of butterfly phy- 
logeny as well as other such studies (de Jong er al., 1996; Caterino et al., 2001) suggested 
that any butterflies of the Papilionoidea could be used as appropriate outgroups to root the 
skipper tree. For this study, three Japanese papilionids, two nymphalids, and two lycaenids 
were collected and used (Table 1). Specimens collected within the past three years and kept 
in 99.596 ethanol or dried (Ochlodes asahinai Shirózu) were used for total DNA extraction. 


DNA extraction 


For most species, total DNA was extracted from an adult individual using SDS-proteinase 
K buffer. Either leg tissue or thoracic muscles were dissected out, squashed in 500 ul pro- 
teinase K solution (10 mM Tris-HCl [pH8.0], 10 mM EDTA, 0.5% SDS, and 0.5 mg/ml 
proteinase K [Boehringer Mannheim GmbH, Germany]), and incubated at 55°C for 2 h. 
The resulting solution was extracted once with phenol saturated with TE buffer (10 mM 
Tris-HCI [pH8.0] and 1 mM EDTA), once with phenol/chloroform/isoamyl alcohol (25: 25: 
1), and once with chloroform/isoamyl alcohol (24: 1). The obtained supernatant was trans- 
ferred to a small centrifuge tube, and RNA was digested with 8 mg/I RNase A (Wako Pure 
Chemical Industries, Ltd, Japan) at 37°C for 30 min. The total DNA was precipitated with 
ethanol, washed with 70% ethanol, dried, and dissolved in 200 ul TE buffer. 


For the remaining species, i. e. Parnara guttata (Bremer & Grey), Thoressa varia (Murray), 
and Erynnis montanus (Bremer), total DNA was extracted with CTAB buffer (Murray & 
Thompson 1980). Whole samples were frozen and ground in liquid nitrogen. The powder 
was put into 600 ul 2xCTAB buffer (2% CTAB, 1.4 M NaCl, 0.1 M Tris-HCl [pH 8.0], 20 
mM EDTA) and incubated at 60°C for 1 h. The solution was then extracted with 700 ul of 
chloroform/isoamylalcohol (25: 1), and the nucleic acid was precipitated with 300 ul 2- 
propanol. The precipitate was rinsed with 300 ul 70% ethanol, dried, dissolved in 50 ul 8 
mg// RNase A, and then incubated at 37°C for 30 min. Total DNA was precipitated with 5 
ul 3 M sodium acetate [pH5.2] and 125 ul ethanol. The precipitate was rinsed with 70% 
ethanol, dried, and dissolved in 200 ul TE buffer. 
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Amplification and sequencing of the ND5 and COI genes 


Total DNA was used as a template for amplification of mitochondrial DNA fragments by 
polymerase chain reaction (PCR). PCR amplification was carried out in 50 ul volume and 
contained 10xEx Taq Buffer, 0.20 mM of each dNTP, 0.8 uM of each primer, 1.3 units Ex 
Taq DNA polymerase (Takara, Japan), and 1 ul DNA solution. The PCR mix was incubated 
at 94°C for 3 min, followed by 35-40 amplification cycles and a final step of synthesis for 7 
min at 72°C. Each amplification cycle included a 1 min period for denaturation at 94°C, a 1 
min period for annealing at 48°C, and a 2 min period for extension at 72°C. Primers used 
for ND5 were V1 (Yagi et al., 1999), KAIL, TIAI, and T2A1. Primers KAIL (GTT 
CTAATATAAGGTATAAATCATAT), TIAL (AATATAAGGTATAAATCATATTA), and T2 
Al (AATATAAGGTATAA ATCATATTAA AQ) were designed by modifications to Al (Yagi 
et al., 1999). Primers used for COI were CO1F1629 (CATTAATTGGAGATGATCAAA 
TTTAT), COIF Ron (GGATCACCTGATATAGCATTCCC), COIF k698 (TACAATTTA 
TCGCCTAAACTTCAGCC), CO1R2952 (ATGTTCAGCAGGAGGAGTATTTTGA), CO 
IR. Pat (TCCAATGCACTA ATCTGCCATATTA), and COIR_Patll (TCCATTACATAT 
AATCTGCCATATTAG). These primer sequences were kindly provided by Dr Ken-Ichi 
Odagiri. The PCR products were separated by electrophoresis on 1.2% SeaKem GTG 
agarose gels (BioWhittaker Molecular Applications, Rockland, ME) and purified with a 
GENE CLEAN III Kit (BIO101) or a Qiaquick Gel Extraction kit, or directly purified with 
a Qiaquick PCR purification kit. Direct sequencing of the purified PCR products was per- 
formed using the Big Dye terminator cycle sequencing kit version 2.0, 3.0, or 3.1 (Applied 
Biosystems, Foster City, CA), with one of the end primers used for amplification, and inter- 
nal primers C2 and A3 (Yagi et al., 1999). The nucleotide sequence was determined with 
an ABI PRISM 377XL DNA Sequencer or ABI PRISM 3100 Genetic Analyzer (Applied 
Biosystems). Nucleotide sequences from a sample were combined with AutoAsssembler 
2.1.1 (Applied Biosystems). 


Phylogenetic analysis 


The nucleotide sequences were aligned using the program fftnsi in MAFFT (Katoh et al., 
2002) version 3.85. 5' and 3' end regions where only a part of the species were determined 
were removed with MacClade 4.06 PPC software (Maddison & Maddison 2000). No sites 
with insertions and/or deletions were found. A 787-bp region in ND5 and 1110 bp of COI 
were subject to phylogenetic analysis. Phylogenetic trees based on nucleotide sequences 
were constructed with the maximum parsimony (MP) method and the maximum likelihood 
(ML) method under the HKY85 model (Hasegawa et al., 1995) by PAUP* version 4.0 beta 
10 (Swofford, 1998). In the MP method, all characters were treated unordered and equally 
weighted. A heuristic search with tree bisection and reconnection (TBR) was performed on 
start trees obtained by random stepwise addition. Random addition was performed 1,000 
times for the original dataset and ten times for the bootstrap datasets in MP analysis, and 
100 times for the original dataset and five times for the bootstrap analysis in ML analysis. 
The confidence level of each branch in the phylogenetic tree was evaluated with the boot- 
strap test (Felsenstein, 1985) based on 10,000 and 100 resamplings for MP and ML analy- 
ses, respectively. Phylogenetic trees based on translated amino acid sequences were per- 
formed with ProtML in the Molphy-2.3b3 package (Adachi & Hasegawa, 1996) under the 
mtRev24 model. Translation was performed with MacClade, with the genetic code set to 
"Drosophila mtDNA." A local rearrangement (nearest neighbor interchange) search was 
performed on the NJ tree obtained with NJDist. For the bootstrap test of the amino acid 
ML tree, 10,000 datasets were prepared with SEQBOOT in PHYLIP 3.573c (Felsenstein, 
1995), and the ML tree was searched for each dataset. The local rearrangement search on 
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some datasets showed oscillation and seemed unlikely to resolve; therefore, only datasets 
that completed within a given time were used for the calculation of bootstrap values with 
CONSENSE. 


Results 


NDS sequences were determined for 37 species of skippers. For the NDS nucleotide se- 
quences used, neither deletions nor insertions were required for multiple sequence align- 
ment. The G+C content (96) was 17.72+0.28 (S.D.) on average and was nearly constant in 
all samples. Since the G+C content strongly deviated from 50%, we used the HKY85 
model (Hasegawa et al., 1995), which does not require an assumption of equal frequencies 
of the four nucleotides, in the likelihood-based analyses. 


Fig. 1 shows the ND5 nucleotide ML tree. In the MP method, of 787 characters, 369 were 
constant and 329 were parsimony informative. A total of four equally parsimonious trees 
with 1,955 steps (CI=0.315, RI=0.428) were obtained. None of the four trees was exactly 
compatible with the conventional classification. Bootstrap values were low in deep branch- 
es, but except for such branches, nucleotide ML, MP, and amino acid ML trees were con- 
gruent. 


In an attempt to obtain stronger support for branches unresolved with NDS gene data, COI 
sequences of 26 skippers and seven outgroup taxa were determined. The CO/ data set alone 
poorly resolved the relationships (data not shown). However, combined analysis of ND5 
and COI genes improved the resolution (Fig. 2) in terms of bootstrap support for some 
groups, although the values were not directly comparable to the analysis of ND5 data alone, 
because the number of taxa was reduced. 


The Coeliadinae were consistently resolved as a monophyletic group in every analysis, 
whereas the Heteropterinae were resolved as monophyletic only in the combined analysis 
(Fig. 2). Tagiades and Pyrgus groups in the Pyrginae were split into two clades. [soteinon 
and other groups in the Hesperiinae were also split into different clades. However, relation- 
ships of these clades and other clades were not resolved with high bootstrap probabilities; 
therefore, monophyly of the Pyrginae or the Hesperiinae was not rejected. 


For genus groups, each Pyrgus group, Tagiades group, Gegenes group, Hesperia group, 
Taractrocera group, and Ancistroides group was consistently resolved as a monophyletic 
group in every analysis. The monophyly of the Gegenes group was only weakly supported 
by the ND5 analysis (52% with nucleotide ML) but was moderately supported by the com- 
bined data set (8296/7796 with nucleotide/amino acid ML). Within the Gegenes group, the 
three species of Parnara formed a clade, but the relationships among Pelopidas, Polytremis, 
and Borbo were unclear. Parnara ogasawarensis was closer to P. guttata than to P. bada, 
with relatively strong support. Within the Hesperia group, Thymelicus diverged first, and 
Hesperia and Ochlodes formed a monophyletic group with moderate to high bootstrap sup- 
port. Within Ochlodes, O. venatus and O. ochraceus formed a monophyletic group support- 
ed by high bootstrap values in the NDS analyses. 


Within the Coeliadinae, Choaspes and Badamia formed a monophyletic group, and Burara 
was basal in the ND5 tree and the ND5+COI tree. The Choaspes-Badamia clade was also 
recovered with COI amino acid analysis (not shown), although it was not found in the COT 
nucleotide analyses, and bootstrap support was low (<50%) in the combined nucleotide 
analyses. 
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52/ 88 99 Pelopidas jansonis Gegenes 
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93/97 Ochlodes venatus 
96/98] 92 Ochlodes ochraceus 
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84/84 Hesperia florinda group 
100/100 r^ Thymelicus sylvaticus 
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Thymelicus lineola 
Suastus gremius Plastingia group 
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Potanthus flavus ] group 
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Aeromachus inachus Astictopterus group 
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92 Udaspes folus group 
52/55 Thoressa varia ーー Astictopterus group 


Isoteinon lamprospilus 
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0.05 substitutions/site 
Fig. l. Maximum likelihood tree of skippers based on nucleotide sequences of ND5 gene (787 bp). 
Bootstrap probabilities by nucleotide ML and MP methods are shown on the branches. 
Bootstrap probabilities by amino acid ML analysis are shown under the branches. The 
branch lengths are proportional to the number of substitutions per site. Higher taxonomic 


categories are indicated to the right. 
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Fig. 2. Maximum likelihood tree of skippers based on nucleotide sequences of ND5 (787 bp) and 


COI (1,110 bp) genes. 


Bootstrap probabilities by nucleotide ML and MP methods are 


shown on the branches. Bootstrap probabilities by amino acid ML analysis are shown under 
the branches. The branch lengths are proportional to the number of substitutions per site. 
Higher taxonomic categories are indicated to the right. 


Discussion 


1) Monophyly of genera 


In this discussion regarding morphological characteristics, especially the male genitalia, we 
refer to the following literature in addition to our own observations: Evans (1937, 1949, 
1951-1955), Shirózu (1960), and Kawazoé & Wakabayashi (1976). 


Our samples included seven genera that were represented by two or three Japanese species. 
These genera were Hasora, Pyrgus, Carterocephalus, Ochlodes, Thymelicus, Pelopidas and 
Parnara. Although each of these genera includes many more species globally, some indica- 
tion of their monophyly may be obtained from our results based on a few Japanese species. 
Bootstrap values of the clades were 100% or nearly 100% for every one of these genera 
(Fig. 1). Therefore, our molecular data strongly support close relationships between or 
among Japanese congeners and their inclusion in each genus as a monophyletic group. 


Considering the support from the molecular data, it is interesting to evaluate some morpho- 
logical characters of the above genera. With regard to the morphology of male genitalia, 
each genus has unique character state(s) as follows: Pyrgus (valval ampulla with a ventro- 
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proximal and a dorsodistal projection), Carterocephalus (elongation of uncus, gnathos, and 
saccus), Ochlodes (branch-like accessory process of aedeagus), Thymelicus (elongation of 
saccus), Pelopidas (a pair of setiferous papilla-like processes on the basal portion of the 
uncus), and Parnara (eaves-like process from the posterior margin of the tegumen, crest- 
like lateral processes on the subapical portion of uncus). These character states may be 
considered specializations from generalized hesperiid genitalia, and are absent or rarely 
found among genera belonging to the same group, except for the complicated aedeagus of 
Ochlodes which is also seen in its allied genera (Burns, 1992). If these character states are 
autapomorphic, then together with molecular data they strongly support monophyly of these 
genera. In Pyrgus, although P. malvae (Linnaeus) is isolated from other congeners because 
of the widely bifurcate uncus of its male genitalia, it was clustered together with P. macula- 
tus (Bremer & Grey) in a clade with 100% bootstrap support in the ND5 gene tree. This re- 
sult strongly supports the monophyly of Pyrgus. However, for Carterocephalus, Ochlodes, 
and Pelopidas, further investigation of molecular data from more congeners (and of allied 
genera of Ochlodes) is necessary to test the monophyly of these genera, because their intra- 
generic morphological diversity is notable, and the samples used in this study do not 
sufficiently cover this diversity. 


2) Intrageneric relationships 


Intrageneric relationships can only be discussed for genera of which at least three Japanese 
species were included in this study, i. e. Ochlodes, Thymelicus, Pelopidas, and Parnara. 


a) Ochlodes. In our analyses, O. venatus (Bremer & Grey) and O. ochraceus (Bremer) 
were more closely related to each other than to O. asahinai, and the bootstrap value of the 
clade consisting of the former two species was extremely high. With regard to morphology, 
sexual dimorphism in wing shape and markings are more strongly developed in O. venatus 
and O. asahinai than in O. ochraceus, and the elongation and curvature of the male genital 
uncus are stronger in the former two. Thus, O. venatus and O. asahinai superficially resem- 
ble each other more than they do O. ochraceus. On the other hand, O. asahinai has its own 
genital specializations, in that the uncus is extremely deeply bifurcate and the gnathos is 
completely fused with the genital dorsum. Ochlodes venatus and O. ochraceus are more 
generalized in these two characters, and O. ochraceus and O. asahinai are more generalized 
in the shape of the valva. However, these two resemblances do not support a close relation- 
ship of each pair. When we simply rely on morphological characters, O. venatus and O. 
asahinai are more closely related to each other than to O. ochraceus. However, the topolo- 
gy of these three species obtained from molecular analyses does not support this specula- 
tion based on morphology. 


b) Thymelicus. Based on our results, T. sylvaticus (Bremer) and T. leoninus (Butler) were 
more closely related to each other than to T. /ineola, and the bootstrap value of the clade 
consisting of the former two species was 100%. Differences in male genitalia and wing 
markings are slight among Thymelicus species. The former two species have similar mark- 
ings on the upper surface of the forewing that are distinctly separated by narrow blackish 
veins. This character state is unique in Thymelicus. If this state is apomorphic, the topolo- 
gy obtained based on wing markings is in accordance with the molecular topology. 


c) Pelopidas. In the molecular analyses, P. mathias (Fabricius) and P. agna (Moore) were 
more closely related to each other than to P. jansonis (Butler), and the bootstrap value of 
the clade consisting of the former two species was 7496 (Fig. 1). Pelopidas mathias and P. 
agna strongly resemble each other in terms of wing markings, including the homonomous 
appearance of postdiscal spots on the hindwings, the linear male sexual brand on the upper- 
side of the forewings, and the complicated cornuti of the male phallus. At least the shape of 
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the sexual brand seems to be apomorphic. On the other hand, P. jansonis has unevenly de- 
veloped postdiscal spots and a large discoidal spot on the underside of its hindwings, and 
the male has a broad inconspicuous sexual brand. The phylogenetic relationship of the 
three species based on the morphological characters agrees well with that determined from 
our molecular data. 


d) Parnara. In the molecular analyses, P. guttata and P. ogasawarensis Matsumura were 
more closely related to each other than to P. bada, and the bootstrap value of the clade con- 
sisting of the former two species was 69%. Parnara ogasawarensis is endemic to the Bonin 
Islands, while the two other species are widely distributed in the Oriental Region, and reach 
Japan. The former was once treated as a subspecies of P. guttata (Evans, 1949), which is 
now known to be distributed in the Bonin Islands, but there are slight differences in the male 
genitalia between the two species, in addition to the conspicuous differences in wing mark- 
ings and body size. The pattern of wing markings in P. ogasawarensis may be interpreted 
as a simple exaggeration and central shift of the P. guttata markings, which are character- 
ized by posteriorly widened postdiscal series in both the fore- and hindwings. Parnara 
bada (Moore) differs from the other two species in the weakly developed wing markings 
and the arched dorsal margin of the male valval ampulla. The phylogenetic relationships in- 
ferred from the morphological data support the molecular topology. 


3) Validity of genus groups proposed by Evans 


Our samples included the following genus groups that contained two or more Japanese gen- 
era in the subfamilies Pyrginae and Hesperiinae (here Evans included the Heteropterinae) 
proposed by Evans (1949): Tagiades gp, Pyrgus gp. Heteropterus gp, Astictopterus gp, 
Ancistroides gp, Plastingia gp, Hesperia gp, Taractrocera gp, and Genenes gp. In the fol- 
lowing section, we discuss the validity of the above-mentioned genus groups based on our 
molecular data. It should be noted that Evans proposed his genus groups based on normal 
taxonomic characters, and at present most of them are not characterized by the cladistic 
method. 


a) Tagiades group. This group is pantropical in its distribution and is represented by nearly 
20 genera, including Tagiades and Daimio, which were analyzed in this study. The boot- 
strap value of the Tagiades-Daimio clade was 10096 (Fig. 1), and clearly supported a close 
relationship between the two genera. However, they are extremely similar to each other in 
many morphological and biological features, and as a consequence our data do not offer any 
support for the Tagiades group as a whole, which includes many diversified genera. 


b) Pyrgus group. This group contains many genera closely related to Pyrgus, and Erynnis is 
very isolated from them. In order to accommodate Erynnis and allied genera in the New 
World, Evans (1951—5) established the Erynnis group. Although the bootstrap value sup- 
porting a sister relationship between the two genera was <50% in the ND5 tree, bootstrap 
support in the combined tree was moderately high, supporting the monophyly of this group. 


c) Heteropterus group. This group includes many genera distributed in the Afrotropical, 
Nearctic, and Neotropical Regions. The two Japanese genera of the group, Leptalina and 
Carterocephalus, did not form their own clade in the ND5 analysis, but their monophyly 
was moderately supported by the combined analysis. Leptalina and Carterocephalus are 
morphologically fairly isolated from each other within the group. Evans (1949) initially 
combined them into one group, but then (1955) created the Carterocephalus group for 
Carterocephalus and the New World genera resembling it. 


d) Astictopterus group. This group includes some Oriental genera in addition to Thoressa 
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and Aeromachus. Evans (1949) classified them into different subgroups: Thoressa in the 
Halpe subgroup and Aeromachus in the Ampittia subgroup. In this study, the two genera 
did not form a single clade but were widely isolated. 


e) Ancistroides group. This group includes more than ten other Oriental genera in addition 
to Udaspes and Notocrypta. The two genera together formed a single clade supported by a 
bootstrap value of 97%. Among the genera in the Ancistroides group, the two investigated 
in this study are extremely similar to each other in terms of male genitalia, morphology of 
immature stages, and larval food plants, although the hindwing markings are quite different 
between the two groups. Moreover, hybrids between U. folus (Cramer) and N. curvifascia 
(C. & R. Felder) have been reported (Yoshida, 2001). Therefore, our results do not offer 
strong support for the monophyly of the whole Ancistroides group. 


f) Plastingia group. This group includes some Oriental genera in addition to Erionota and 
Suastus. Evans (1949) classified Suastus in the Plastingia subgroup and Erionota in the 
Erionota subgroup. In our analysis, the two genera were included in different clades, for 
which bootstrap values were «5096. The molecular data thus do not offer any support for 
the monophyly of this group. 


g) Hesperia group. This group contains some Holarctic genera, most of which strongly re- 
semble each other. Our samples included Hesperia, Ochlodes, and Thymelicus, which 
formed a single clade with a bootstrap value of 8496. Hesperia and Ochlodes resemble 
each other, but Thymelicus is morphologically rather isolated from the former two genera; 
however, the bootstrap values support the monophyly of the entire Hesperia group. In 
terms of genital morphology, this group is characterized by the deeply incised anterior mar- 
gin of the tegumen, which may be considered autapomorphic to the group. 


In this study, Ochlodes and Hesperia together formed a clade, which was the sister group of 
Thymelicus. Thymelicus is unique in having an elongate specialized saccus and a single 
generalized uncus in the male genitalia, whereas in Ochlodes and Hesperia, the saccus is 
normally short but the uncus is bifid distally. The latter character state seems to be synapo- 
morphic, combining the two genera. Therefore, the tree topology of the three genera is 
identical whether based on the morphological or molecular data. 


h) Taractrocera group. This group includes many Oriental genera which superficially re- 
semble each other. It is morphologically characterized by the undivided apical portion of 
the valva and the absence of a male genital gnathos, supporting the monophyly of this 
group. In our analysis, Telicota and Potanthus together formed a clade with a bootstrap 
value of «50926; however, the combined tree supported the monophyly of the group, with 
97% bootstrap support for amino acid sequences. 


1) Gegenes group. This group includes many temperate to tropical genera in the Old World, 
which are similar to each other in the dark brown upperside and yellowish brown underside 
of the wings with subhyaline whitish spots, as well as in the strongly swollen antennal club 
with a short apiculus. This group was represented in our samples by Polytremis, Pelopidas, 
Borbo and Parnara. Although these genera formed a clade, its bootstrap value was only 
5896. This value suggests diversity in the group, as the structure of the male genital dorsum 
is extremely diversified among the included genera, despite the above-mentioned external 
resemblances. In our analysis, Pelopidas and Borbo formed a clade supported by a low 
bootstrap value; the sister group of this clade was Polytremis. The clade formed by these 
three genera was supported by a high bootstrap value, and it is morphologically supported 
by the bifurcate uncus of the male genitalia. Parnara differs considerably from the three 
genera in the structure of the genital dorsum, which has a teguminal eave and a single apex 
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of the uncus with a pair of crest-like processes on it. The topology obtained by the molecu- 
lar data suggests that further morphological analysis of these intergeneric relationships is 
required. 


4) Classifications higher than the genus group level 


Because our samples were confined to Japanese species and the NDS gene is rather short, 
we did not aim to test or evaluate hesperiid classification higher than the genus group. 
However, our samples covered four subfamilies, and the resulting trees were in some degree 
related to higher categories; therefore, a few brief comments are warranted. 


a) Intergroup relationships of the Hesperiinae. In his work, Evans did not comment on any 
strict phylogenetic relationships among his genus groups in the subfamily Hesperiinae. In 
our analysis, the bootstrap values of all branches combining different genus groups were 
«8090, and no information was obtained about the relationships among the genus groups 
from our analysis of the NDS gene. 


b) Monophyly of subfamilies. Topology and bootstrap values differed by subfamily, as dis- 
cussed below. 


i) Coeliadinae. In both MP and ML trees, the genera of the subfamily Coeliadinae formed a 
single clade, and the bootstrap value was moderately high in the combined analysis, sup- 
porting monophyly of this subfamily. This subfamily is morphologically well characterized 
by some unique (probably apomorphic) character states, including the very slender style- 
like third labial palpomere, which is projected forward and densely covered with appressed 
scales. Intergeneric relationships in the subfamily inferred from molecular data are also in- 
formative. The close relationship between Choaspes and Badamia was supported with a 
bootstrap value of 86% in the ML analysis using combined amino acid data, whereas boot- 
strap support from the nucleotide data was <50%. The topology obtained from the ML tree 
showed a close relationship between Hasora, Choaspes, and Badamia, with a 54% boot- 
strap value. 


ii) Pyrginae. In this study, two groups, Pyrgus gp and Tagiades gp, were included in the 
subfamily Pyrginae. As already discussed, monophyly is well supported in our analyses for 
Pyrgus gp but poorly supported for Tagiades gp. However, the subfamily became para- 
phyletic in the combined tree, as the Tagiades group was clustered with the subfamily 
Coeliadinae. Consequently our analyses do not offer any valid information on the systemat- 
ic position and monophyly of the subfamily Pyrginae. 


iii) Heteropterinae. In this study, two genera, Leptalina and Carterocephalus, were includ- 
ed in the subfamily Heteropterinae. In the combined tree, the two genera formed their own 
cluster, but its bootstrap value was <70%. Therefore, our analysis is not fully informative 
as to the potential monophyly and systematic position of the subfamily Heteropterinae. 


iv) Hesperiinae. Most of the hesperiine genus groups, except for Thoressa and Isoteinon, 
formed a single large cluster, although its bootstrap value was <60% in the combined tree. 
In the ML tree, the cluster included Carterocephalus (Heteropterinae), and the /soteinon- 
Thoressa branch was isolated from all other genera investigated in this work. Therefore, 
our analysis does not offer valid information on the monophyly of the subfamily 
Hesperiinae. 
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摘 要 


ミト コン ドリ ア APD? お よび CO7 遺 伝 子 の 塩基 配列 に よる 日 本 産 セ セリ チョ ウ の 分 子 糸 統 
(百々 由 項 子 ・ 三 枝 豊 平 ・ 千 葉 秀 辛 ・ 西 山 智明 ・ 広 渡 俊哉 ・ 石 井 実 ・ 八 木 其 司 ・ 長 谷 部 光 
泰 ・ 毛 利 秀 雄 ) 


日 本 産 セ セリ チョ ウ 全 37 種 に つい て , ミト コン ドリ De F CO/ 遺伝 子 領域 
1.110 塩基 対 を 用 いて , PSHM (MP) on: (ML) に より 分 子 系 統 樹 を 構築 し た . 本邦 に 2 種 以 上 を 
産 する 7 属 ,. す な わ ち ビ ロ ウ ド セ セ リ 属 , チ ャ マダ ラ セ セ リ 属 , タ カネ キマ ダラ セセリ 属 , コキ マダ ラ セ 
セリ 属 , スジ グロ チャ バネ セセリ 属 , チ ャ バネ セセリ 属 お よび イチ モン ジ セ セリ 属 に つい て は . 
て ブー トス トラ ッ プ 値 が 100 か それ に 近く , 単 系 統 性 が 支持 され た . 本邦 に 3 種 以 上 を 産 する 4 属 の 
内 の 種 問 関係 は 次 の 通り 0. コ キマ ダラ セセリ 属 で は , E ec p UE 
と アサ ヒナ キマ ダラ セセリ より 近 縁 . スジ グロ チャ バネ セセリ 属 で は , ス ジグ ロ チ ャ バネ セセリ と ヘリ 
グロ チャ バネ セセリ が それ ら と カラ フト セセリ より 近 縁 . チ ャ バネ セセリ 属 で は , チ ャ パネ セセリ と ト 
ガリ チャ パネ セセリ が それ ら と ミヤ マチ ャ バネ セセリ より 近 縁 . イチ モン ジ セ セリ 属 で は , イチ モン ジ 
セセリ と オガ サワ ラ セ セ リ が それ ら と ヒメ イチ モン ジ セ セリ より 近 緑 . 外部 形態 に よっ て 推定 され た 
Evans の 属 グ ルー プ の 系 統 推定 は 概ね 支持 され た . 
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